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ABSTRACT 

A  summary  of  Bragg  deflection  theory  and  various  approaches 
to  direct  detection  acoustooptic  spectrum  analysis  (AOSA) 
modeling  is  presented.   A  suitable  model  is  chosen  and 
extended  to  include  the  effects  of  diffraction  efficiency, 
transducer  efficiency,  irradiance  profiles  of  incident  laser 
illumination,  aperture  size  of  the  Bragg  cell,  and  the  acoustic 
attenuation  experienced  by  the  acoustic  wavetrain  generated  by 
the  input  r-f  signal.   A  Fortran  program  is  developed  to  model 
the  AOSA  and  predict  the  output  image  plane  intensity  profiles. 
A  second  version  of  the  program  includes  a  time  variable  per- 
mitting dynamic  simulation  of  the  system  response. 
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I.   INTRODUCTION 

A.   OVERVIEW 

Assessment  of  emergent  technologies  leads  to  a  continuing 
focus  on  acoustooptical  techniques  for  spectral  analysis  and 
filtering/excision  [Ref .  1] .   Integrated  optics  with  distinc- 
tive advantages  in  size  and  weight  easily  meet  the  require- 
ments for  radar  warning  receivers.   Bulk  wave  devices  offer 
superior  dynamic  range,  bandwidth,  and  resolution  over  inte- 
grated optical  devices.   Foreseeable  applications  include 
surveillance  and  electronic  intelligence. 

An  acoustooptical  spectrum  analyzer  transforms  a  wideband 
time  domain  signal  into  an  optical  spatially  detected  fre- 
quency domain  signal.   Removal  of  narrowband  interference 
may  be  accomplished  by  optical  excision  techniques,  to  selec- 
tively notch  filter  unwanted  spectral  components.   U.  S. 
military  applications  in  this  area  could  greatly  benefit 
from  the  speed  of  response  and  high  quality  interference  re- 
jection capabilities  possible  with  these  techniques.   The 
potential  of  these  techniques  was  recognized  and  led  to  the 
establishment  of  a  joint  NAVELEX  350/DARPA  program  to  examine 
these  areas. 

The  technical  objectives  of  this  effort  were: 

1.  Survey  potential  applications  and  analyze  technical 
requirements  imposed  by  each. 

2.  Analyze  alternative  optical  excision  configurations 

and  components. 
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3.  Develop  and  procure  major  system  components. 

4.  Fabricate  and  measure  alternative  broadband  excisor 
configurations  to  verify  analytical  results. 

5.  Fabricate  and  field  test  a  brassboard  model. 

6.  Design  a  preliminary  advanced  development  model. 
The  principal  civilian  contractor  for  this  project  is  PROBE 
Systems,  Inc.   PROBE,  assisted  by  the  developmental  efforts 
at  the  Naval  Research  Laboratory  (NRL)  and  Naval  Postgraduate 
School  (NPS) ,  was  tasked  with  the  attainment  of  these  objec- 
tives.  NRL's  specific  program  assistance  was  the  development 
of  a  photodichroic  optical  clipper  to  be  evaluated  by  PROBE. 
The  research  summarized  in  this  thesis  is  a  continuation  of 
research  initiated  as  part  of  the  NPS  contribution  to  the 
development  program. 

B.   THESIS  OBJECTIVE 

The  objective  of  this  work  effort  was  the  attainment  of 
an  understanding  of  the  principles  of  direct  detection  acousto- 
optical  spectrum  analysis .   Previous  work  at  the  Naval 
Postgraduate  School  resulted  in  the  construction  by  Regan 
[Ref.  2]  of  a  working  laboratory  acoustooptical  spectrum 
analyzer  and  the  subsequent  modification  of  the  existing 
equipment  by  Smith  [Ref.  3]  to  function  as  an  acoustooptical 
spectrum  excisor  (AOSE) .   The  acoustooptical  spectral  excisor 
utilized  optical  excision,  the  effect  of  an  object  placed 
in  that  part  of  the  optical  path  where  the  undesirable 
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frequency  components  have  been  spatially  spread.   The  optical 
excision  research  included  propagation  considerations  and 
diffraction  by  the  obstacle.   The  effectiveness  of  the  optical 
excisor  in  removing  interference  is  dependent  upon  these 
factors  and  more  fundamentally  upon  the  ability  of  the  de- 
flected beam  irradiance  pattern  to  accurately  represent  the 
input  r-f  signal  spectral  components. 

C.   ACOUSTOOPTICAL  SPECTRAL  ANALYSIS 

An  acoustooptical  spectrum  analyzer  (AOSA)  consists  of 
an  acoustooptical  light  deflector  (a  Bragg  cell) ,  a  coherent 
light  source  (a  laser) ,  and  a  detector.   An  input  r-f  signal 
excites  a  piezo-electric  transducer  mounted  on  the  acoustic 
medium  on  the  Bragg  cell.   This  launches  an  acoustic  wave- 
train  the  length  of  the  cell  which  modulates  the  refractive 
index  of  the  interaction  medium  by  the  elasto-acoustic 
process.   The  modulated  interaction  material  functions  as 
a  moving  diffraction  grating  for  the  incident  laser  light. 
A  one-dimensional  angular  pattern  of  diffracted  beams  is 
produced  and  focused  on  a  photodetector  array.   The  pattern 
replicates  the  spectral  content  of  the  original  r-f  signal. 
The  intensity  of  the  individual  beams  is  directly  proportional 
to  the  r-f  signal  component  intensity.   The  deflection 
angles  of  the  individual  beams  are  proportional  to  the  r-f 
signal  frequency  component.   The  detector  array  photosites 
video  detect  the  incident  beam  patterns  and,  in  a  square  law 
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fashion,  convert  the  spectral  information  for  electronic 
transmission  to  a  display  oscilloscope.   Figure  1  shows  the 
NPS  laboratory  AOSA.   Figure  2  shows  a  typical  oscilloscope 
display  output  from  the  system.   Korpel  [Ref.  4]  includes  a 
basics  discussion  on  the  theory  of  Bragg  cell  applications 
for  acoustooptical  spectrum  analysis  in  an  overview  on 
acoustooptics . 

D.   SPECIFIC  GOALS 

There  were  three  initial  goals  for  this  work.   The  first 
was  the  study  of  the  factors  affecting  the  Bragg  cell 
deflection  of  incident  laser  light.   The  second  was  the 
analysis  of  the  effect  of  these  factors  on  the  replication 
of  the  input  r-f  signal  spectrum  in  the  deflected  beam 
pattern.   The  third  was  the  development  of  a  Fortran  program 
that  would  model  the  acoustooptic  spectrum  analysis  process. 

1.   Bragg  Diffraction 

The  study  and  analysis  of  the  factors  affecting 
Bragg  deflection  was  focused  towards  an  understanding  of  the 
ability  of  the  deflected  beam  pattern  to  faithfully  replicate 
the  spectral  components  of  the  input  r-f  signal.   This 
knowledge  was  considered  essential  in  order  to  construct  a 
viable  computer  simulation  model.   A  study  of  Bragg  diffrac- 
tion and  the  various  approaches  to  modeling  of  the  deflected 
beam  irradiance  profiles  was  undertaken.   Previous  works  at 
NPS  by  Regan  [Ref.  2]  and  Smith  [Ref.  3]  were  used  as 
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Figure  2.   Oscilloscope  Display  of  a  70  MHz  Monotonic  Input 
Signal  to  the  NPS  Acoustooptic  Spectrum  Analyzer 
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starting  points.   Bragg  diffraction  is  traditionally  theo- 
rized by  use  of  complex  coupled  mode  wave  equations.   The 
complexity  of  these  equations  rendered  them  unsuitable  for 
the  first  order  simulation  sought.   A  study  of  other 
approaches  was  conducted  to  arrive  at  a  method  that  would 
provide  both  adequate  accuracy  and  a  reasonable  level  of 
complexity. 

2.   Computer  Model 

An  approach  to  the  modeling  of  the  AOSA  was  selected 
based  on  the  work  of  Hecht  [Ref .  5] .   He  shows  the  image 
plane  intensity  as  the  squared  magnitude  of  the  convolution 
of  the  spatial  Fourier  transform  of  an  illumination  function 
and  the  transform  of  the  converted  acoustic  wave.   Some 
modifications  in  the  form  of  transmission  factors  were  applied 
by  algebraic  multiplication.   Hecht 's  illumination  factor  is 
a  spatial  weighting  function  composed  of  the  optical  ampli- 
tude profile,  the  optical  aperture  size  of  the  Bragg  cell, 
and  the  acoustic  attenuation  of  the  interaction  material. 
The  added  transmission  factors  were  the  diffraction  efficiency 
and  the  transducer  efficiency. 

A  Fortran  program  was  developed  to  model  the  AOSA 
given  the  specific  physical  parameters  of  the  system  and 
desired  simulation  requirements  such  as  laser  wavelength  and 
laser  beam  profile.   This  program  was  used  to  study  the 
quality  of  the  replication  of  the  r-f  signal  spectral  components 
The  major  thrust  of  this  effort  was  to  produce  this  model  for 
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utilization  with  previously  developed  optical  excision 
programs.   While  the  program  is  configured  to  simulate 
various  component  deflectors,  it  was  structured  about  the 
AOD-7  0  model  Bragg  deflector,  manufactured  by  the 
IntraAction  Corp. ,  which  is  scheduled  for  future  installa- 
tion in  the  NPS  laboratory  AOSA.   The  AOD-70  Bragg  deflector 
is  constructed  with  high  quality  flint  glass  and  Indium 
bonded  Lithium  Niobate  stepped  array  transducers.   It  has 
a  50  MHz  to  90  MHz  bandwidth  and  a  center  frequency  of  70  MHz 

During  the  development  of  the  program  it  became 
apparent  that  the  addition  of  a  third  dimension,  time,  to 
the  model  would  significantly  enhance  the  ability  of  the 
model  to  demonstrate  the  response  of  the  AOSA  under  various 
simulation  conditions.   A  second  Fortran  program  was  written, 
essentially  a  variation  of  the  first,  to  provide  this  dynamic 
capability. 

E.   PERSPECTIVE 

The  application  of  both  AOSA  and  AOSE  techniques  in  wide- 
band, high  density  signal  environments  is  dependent  on  the 
accuracy  of  the  diffracted  beam  profile  in  its  replication 
of  the  input  signal  spectra.   The  classical  model  for 
determination  of  the  diffracted  beam  profile  is  presented 
by  Hecht  [Ref .  5] .   In  a  simplified  form  the  diffracted  beam 
profile  is  the  convolution  of  the  spatial  Fourier  transform 
of  the  illumination  field  incident  on  the  Bragg  cell  and  the 
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spatial  transform  of  the  generated  acoustic  wave  within  the 
cell.   In  detail  both  the  illumination  function  and  the 
acoustic  wave  are  subject  to  modifying  factors.   The  signifi- 
cant factors  include  diffraction  efficiency,  transducer 
efficiency,  the  irradiance  profile  of  the  incident  laser 
light,  the  aperture  size  of  the  Bragg  cell,  and  the  acoustic 
attenuation  experienced  by  the  acoustic  wave  as  it  transits 
the  cell.   This  thesis  is  an  attempt  to  provide  a  more 
detailed  model  by  incorporating  these  factors  in  a  first 
order  simulation  to  aid  in  the  study  of  acoustooptic  spectrum 
and  excision  techniques. 
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II.   ACOUSTOQPTICAL  SPECTRAL  ANALYSIS  MODELING 

A.   ACOUSTOOPTICAL  INTERACTION  THEORY 

As  an  acoustic  wave  propagates  in  an  isotropic  material, 
it  produces  periodic  modulation  of  the  index  of  refraction 
by  means  of  the  elasto-acoustic  effect.   The  interaction  of 
incident  light  waves  with  acoustic  waves  is  termed  acousto- 
optical  coupling  and  is  traditionally  modeled  utilizing 
coupled  wave  equations.   The  basic  theory  of  acoustooptical 
interaction  in  isotropic  media  was  well  developed  prior  to 
the  advent  of  the  laser.   An  excellent  discussion  of  early 
theoretical  work  can  be  found  in  Ref.  6  by  Born  and  Wolf. 
The  advent  of  the  laser  stimulated  research  in  acoustooptic 
interactions  as  devices  were  developed  for  the  deflection 
and  modulation  of  incident  laser  light.   Development  of 
superior  acoustooptic  materials  and  highly  efficient  broad- 
band transducers  have  given  impetus  to  rapid  advances  in 
these  acoustooptic  applications. 

Various  models  describing  acoustooptic  interaction  have 
been  proposed.   A  thin  grating  model  (as  it  assumes  a  negli- 
gible depth  for  the  sound  field)  was  proposed  by  Rosenthal 
[Ref.  7],  Korpel  and  associates  [Ref.  8],  and  Lambert  [Ref.  9] 
The  thin  grating  model  is  essentially  a  first  order  solution 
of  the  diffracted  light  problem.   A  second  approach  known  as 
the  physical  optics  model  incorporates  the  Green's  function 
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integral,  as  proposed  by  Gordon  [Ref.  10].   An  extremely- 
rigorous  model,  the  mathematics  involved  are  quite  complex. 
A  third  approach  is  termed  the  phenomenological  geometrical 
optics  model  as  proposed  by  Johnson  [Ref.  11]  and  ranks  midway 
in  complexity  between  the  thin  grating  and  physical  optics 
models. 

For  the  purposes  of  this  work  a  highly  satisfactory 
approach  to  acoustooptical  diffraction  is  found  in  a  series 
of  works  by  Chu,  Kong,  and  Tamir  [Refs.  12-16].   They  have 
considered  diffraction  of  arbitrary  profiled  optical  beams 
by  periodically  modulated  layers  and  have  achieved  closed 
form  solutions  for  the  resulting  diffracted  far-field 
intensity  patterns.   A  summary  of  their  results  is  presented 
below. 

As  shown  in  Fig.  3,  a  bounded  beam  is  incident  on  a 
periodically  modulated  dielectric  layer  with  permittivity  of 
the  form 

e(z)  =  £2(1  +  M  cos  27rz/d)  (1) 

where 

e2  =  the  relative  permittivity  of  the  slab  in  the 

absence  of  modulation 
M  =  the  modulation  index 
d  =  periodicity 
The  slab  has  thickness  L  and  is  bounded  by  a  dielectric 
medium  with  relative  permittivity  £i  for  x^O  and  £3  for  x-L. 
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Figure  3.   Bragg  Diffraction  Geometry  (after  Ref.  12) 
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The  electric  field  of  an  incident  beam  may  be  represented  as 


E.       (x 
mc 


,z)    =.  /  G(SQ)    exp  [ J (?q1]x  +    3Q    z)]   d6Q 


(2) 


with 


3     =    (2tt/X)  VeT  sin   9  (3) 


^x)    =    (2TT/X)VrT  cos 


-J  (2tt/\)2    ei-    602  (4) 

G($  )  is  the  angular  spectral  amplitude  of  the  incident  beam 
at  the  entrance  plane  x=0. 


G(6o>  ~h     J   Einc  (0'z)  exP  M  V  dz 

-  00 

=  7{Einc  (0'z)}  (5) 

Letting 

E.    (0,z)  =  F(z)  exp  (jbz)  (6) 

with 

(1) 


b  =  (2tt/X)  ve  i  sin 

=  (2tt/A)  \/Tl    sin   9^3)  (7) 

where 

6    =  angle  of  incidence  of  beam  axis 
o       3 

(3) 
0    =  angle  of  refracted  beam 
o 

F(z)  =  beam-profile  at  x=0 

then  Equation  (5)  becomes 

G(8o)  =  j^    J   F(z)  exp  [-j (3o-b)zJdz  (8) 
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Various  beam  profiles  F(z)  and  their  corresponding  angular 
spectral  functions  G(3  )  are  listed  in  Table  I. 

For  an  optical  beam  at  exactly  Bragg  incidence  we  have 
e(1)  =  e'1'  and  e<3>  =  e  <3»  ,  and 

O  B  O  B 


,(1) 

B 


b  =    (2tt/X)V£T     sin 
=    (2TT/X)\/rr     sin 

=  7T/d 


,(3) 

B 


The  far  field  pattern  for  the  Bragg  deflected  beam  is  given 

by 

P(6)  =  G(<|>)  T(<j>)  (10) 

where 

<t>  ■  (2tt/X)  VTT  sin  6  +  ir/d 

T(<£)  =  is  the  complex  transmission  coefficient. 
A  description  of  this  coefficient  may 
be  found  in  Appendix  C. 

Thus  the  far-field  pattern  for  a  Bragg  deflected  beam  is 

essentially  the  Fourier  transform  of  the  aperture  field  at 

the  boundary  of  the  layer  modified  slightly  by  a  transmission 

factor. 

While  the  results  obtained  by  Chu,  Kong,  and  Tamir  are 

extremely  useful  in  approaching  the  construction  of  a  model 

for  the  direct  detection  acoustooptical  spectrum  analyzer/ 

they  are  predicated  on  the  periodic  modulation  of  the  layer 

by  a  single  frequency  acoustic  wave. 
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B.   ACOUSTOOPTIC  SPECTRUM  ANALYSIS  THEORY 

The  acoustooptic  spectrum  analyzer  consists  of  only  four 
components;  a  light  source,  Bragg  cell,  transform  lens,  and 
a  photodetector  array.   The  basis  for  acoustooptic  spectrum 
analysis  is  the  frequency  dispersion  of  acoustooptic  diffrac- 
tion.  Essentially  each  frequency  component  of  the  input 
signal  will  cause  a  deflection  of  the  incident  light  beam. 
A  typical  interaction  geometry  is  shown  in  Fig.  4  [Ref.  17]. 
A  summary  of  the  description  of  this  interaction  geometry 
from  Ref.  17  is  presented. 

Two  discrete  frequencies  are  shown  with  corresponding 
acoustic  wavelengths  A  i  and  A2.   These  propagate  in  the  x 
direction  with  velocity  V  .   The  medium  is  optically  iso- 
tropic with  an  index  of  refraction  n.   The  thickness  of  the 
ultrasonic  grating  along  the  x  axis  is  L.   A  collimated 
beam  of  light  of  freespace  wavelength  A  is  incident  on  the 
sound  field  at  an  angle  9.  from  the  z  axis  in  the  xy  plane. 
A  sound  wave  at  frequency  f  interacting  with  the  incident 
light  wave  will  generate  a  principal  diffracted  beam  separated 
from  the  incident  beam  by  an  angle  8 .  +  9 ,  by  the  well  known 
Bragg  law: 

-1  /  Xf 


+  8,  =  2  sin 


id  \  2nV 


nVs 


(11) 
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Figure  4.   Bragg  Cell  Operation 
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where  8.  +  9,  <      0.1  rad 
l    d 

External  to  the  medium 

Thus  the  diffraction  angle  is  approximately  proportional  to 
the  input  frequency.   Figure  5  shows  a  classical  optical  ray 
trace  diagram  for  the  acoustooptical  spectrum  analyzer  con- 
figuration.  Figure  6  shows  the  corresponding  optical  plane 
schematic.   In  the  output  transform  plane,  position  is 
approximately  related  to  frequency  by 

X'^)F  (13) 

where 

F  =  transform  lens  focal  length. 

At  this  point,  by  considering  the  input  signal  as  a 
composite  of  frequencies  each  of  which  will  cause  a  deflected 
beam  in  accordance  with  the  presented  interaction  theory, 
an  approach  to  the  construction  of  a  model  may  be  undertaken. 

Schulman  [Ref.  18]  and  Vatz  [Ref.  19]  have  proposed  a 
set  of  optical  block  diagram  symbols  which  show  the  process 
elements  and  indicate  the  mathematical  function  each  element 
performs.   Figure  7  shows  the  AOSA  using  these  process 
elements. 

The  laser  light  source  SL,  and  beam  expander  lenses,  T„ 
and  T  ,  are  shown  as  providing  an  input  to  the  Bragg  cell 
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Figure  5.   Classical  Ray  Trace  Diagram  for  an  Acoustooptical 
Spectrum  Analyzer 
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Figure  6.   Optical  Plane  Schematic  for  an  Acoustooptical 
Spectrum  Analyzer 
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with  transverse  light  distribution  a(x)  of  general  form,  an 
arbitrary  profile  modified  by  spatial  noise,  lens  aberrations, 
and  Bragg  cell  aperture.   The  input  signal  f(t)  is  converted 
to  an  acoustic  wave  in  the  E/0  block  which  represents  the 
piezo-electric  transducer.   This  acoustic  traveling  wave 
f(x-vt)  propagates  through  the  interaction  region. 

In  the  Bragg  cell  interaction  the  light  modulation  can 
be  considered  a  multiplication  process  in  which  the  diffracted 
light  power  is  proportional  to  the  product  of  the  light  inten- 
sity and  the  acoustic  wave  power. 

i(x)  =  a(x)  f(x-vt)  (14) 

where 

x  =  v  t  (15) 

The  transform  lens  T  produces  the  Fourier  transform. 

I(s)  =^{a(x)  f  (x-vt)  } 

=  A(s)  *  F(s)  exp  [j2tt  S(vst)]  (16) 

The  transform  is  the  convolution  of  the  Fourier  transform 
of  the  illumination  function  and  the  transform  of  the  acoustic 
wave.   The  exponential  factor  is  a  Doppler  frequency  shift. 
As  the  detector  is  a  square  law  device,  the  output  from  the 
AOSA  is  the  square  of  the  magnitude  of  the  intensity  function. 

|A(s)  *  F(s)  |2  (17) 

It  has  been  shown  that  the  output  of  the  AOSA  is  the 
square  of  the  magnitude  of  the  intensity  function.   The  in- 
tensity function  may  be  considered  as  the  convolution  of  the 
spatially  transformed  illumination  function  and  the  spatially 
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transformed  acoustic  wave.   Hecht  [Ref.  20]  treats  the  illumi- 
nation function  as  a  spatial  weighting  function  composed  of 
a  number  of  factors  including  optical  amplitude  profile, 
aperture  size,  and  acoustic  attenuation.   The  AOSA  frequency 
plane  intensity  distribution  would  thus  be  proportional  to 
the  squared  amplitude  of  the  Fourier  transform  of  the  weighting 
function  with  respect  to  the  spatial  frequency. 


/oo 
W(U 


F(s)  =   /   W(u)  exp  (-j2?rsu)du  (18) 


where  u  =  [—  -  •*•] 

T     2 


t  =  finite  period  of  signal  analyzed 
S  —  f T  -  the  normalized  frequency 

Further  refinement  in  the  basic  model  is  made  by  the 

consideration  of  the  effects  of  transducer  efficiency  and 

diffraction  efficiency.   These  factors  modify  the  frequency 

components  in  the  acoustic  wave. 

C.   AOSA  MODEL  CONSTRUCTION 
1.   Irradiance  Profiles 

The  choice  of  the  optical  beam  profile  incident  on 
the  Bragg  deflector  is  a  complex  decision.   It  would  appear 
desirable  to  fill  the  entire  deflector  aperture  with  a  uniform 
intensity  light  below  the  saturation  level  of  the  AOSA.   The 
irradiance  profile  in  this  case  would  be  that  of  a  plane  wave. 
A  uniform  intensity  filling  the  entire  aperture  will 
maximize  acoustooptic  interaction  and  provides  maximum 
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incident  light  to  be  contributed  to  the  deflected  beam 
components.   Abrupt  windowing  will  truncate  the  beam  profile 
and  cause  sidelobes  to  appear.   These  sidelobes,  if  they  are 
of  sufficient  magnitude,  may  cause  erroneous  components  to 
appear  at  the  image  plane  or  mask,  low  intensity  components. 

To  reduce  the  sidelobe  levels,  a  less  discontinuous 
transition  in  the  absorption-transmission  regions  is  necessary 
A  suitable  choice  would  be  a  Gaussian  irradiance  profile  with 
a  spherical  phase  wavefront.   This  irradiance  profile  common- 
ly occurs  as  a  laser  output  and  has  well  understood  propaga- 
tion behavior.   Smith  includes  a  handy  propagation  program  as 
part  of  his  work  [Ref .  3] .    This  program  calculates  Gaussian 
beam  phase  front  radius  and  spot  size  for  propagation  through 
any  number  of  thin  lenses  spaced  at  arbitrary  distances. 
Since  the  Gaussian  incident  beam  profile  would  not  fill  the 
entire  aperture  with  uniform  intensity  and  assuming  continued 
operation  below  the  AOSA  saturation  level,  less  light  inten- 
sity would  be  available  for  diffraction  into  representative 
beams.   This  will  result  in  the  loss  of  sensitivity  and  the 
possibility  of  an  inability  to  detect  low  amplitude  frequency 
components  in  the  presence  of  strong  frequency  components . 

Other  beam  irradiance  profiles  are  possible  and 
might  be  employed  for  their  specific  deflected  beam  charac- 
teristics.  Figure  8  taken  from  Ref.  15  shows  the  deflected 
beam  profiles  for  various  incident  profiles.   Such  curves 
are  useful  for  comparison  of  expected   sidelobe  levels. 
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Figure  8.   Bragg  Deflection  Beams  for  Various  Incident 
Irradiance  Profiles  (after  Ref .  15) . 
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Hecht's  [Ref.  5]  implementation  of  the  Gaussian 
irradiance  profile  is  of  the  form: 


wi (t,T,T)  =  exp 


-  M^-\ 


(19) 


where 


T   =  D/2W 

D  =  effective  aperture 

W  =  Gaussian  waist  size 
o 

A  significant  design  question  in  the  utilization  of  Gaussian 

irradiance  profiles  is  the  choice  of  the  ratio  of  the  aper- 

-2  . 

ture  to  the  e   intensity.   This  is  discussed  also  by  Smith 

[Ref.  3] .   Beam  spot  size  is  defined  as  the  perpendicular 

distance  from  the  beam  centerline  to  the  point  where  the 

-2 

irradiance  has  dropped  by  a  factor  of  e   .   To  use  the  full 

information  capacity  of  the  Bragg  cell,  the  spot  size  must 
be  enlarged  to  illuminate  the  complete  cell  aperture.   In 
this  situation  the  aperture  will  severly  truncate  the  Gaussian 
beam  resulting  in  a  profile  similar  to  that  of  a  sine 
function.   Varying  the  ratio  of  the  spot  size  to  the  aperture 
will  alter  the  level  of  the  sidelobes  of  the  profile.   Smith 
limited  his  spot  size  to  avoid  sidelobes  and  thus  isolate 
their  effect  in  his  study  of  optical  excision. 

Hecht  expresses  his  spatial  weighting  functions  in 
equivalent  time  variables.   The  linear  relationship  between 
the  spatial  and  time  variables  and  the  analogous  relationship 
between  the  spatial  frequency  and  frequency  provide  great 
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flexibility  in  the  construction  of  the  computer  program  to 
model  the  AOSA.   The  Fourier  and  inverse  Fourier  transforms 
utilized  in  the  AOSA  program  make  no  distinction  as  to  the 
transforming  variables. 

Only  plane  and  Gaussian  irradiance  profiles  were 
included  in  the  AOSA  program.   It  would  be  a  simple  matter 
to  incorporate  other  profiles  such  as  the  Lorentzian  or 
triangle  beams.   However,  the  frequency  with  which  the  more 
exotic  profiles  would  be  encountered  does  not  merit  the 
addition.   Figure  9  shows  the  comparative  image  plane 
profiles  for  both  irradiance  profiles  with  a  monotonic  70  MHz 
input  signal.   In  Figure  9  the  plane  wave  illumination  chosen 
completely  illuminated  the  aperture  of  the  Bragg  cell.   The 
Gaussian  illumination  used  an  arbitrarily  chosen  3.8  milli- 
meter spot  size.   Comparison  shows  a  near  50%  reduction  in 
the  level  of  the  Gaussian  sidelobes  to  the  level  of  the  plane 
wave  sidelobes  but  with  a  tradeoff  of  a  decrease  in  the  peak 
intensity  response.   The  units  for  the  Image  Plane  Intensity 
axis  in  this  figure  and  all  succeeding  AOSA  program  figures 
are  arbitrary. 

Figure  10  shows  the  image  plane  intensity  for  an 
incident  Gaussian  beam  with  a  spot  size  of  1.945  millimeters. 
This  is  the  "optimum"  spot  size  according  to  Smith's 
criteria  for  sidelobe  suppression.   The  sidelobes  are  barely 
visible  but  again  a  tradeoff  results  in  an  additional  decrease 
in  peak  intensity.   Use  of  a  Guassian  beam  with  this  spot 
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Produced  by  a  Monotonic  7  0  MHz  Signal  with  an 
Incident  Gaussian  Beam  and  a  Simple  Transducer 
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size  narrows  the  portion  of  the  Bragg  cell  aperture  receiving 
full  illumination.   The  portion  receiving  the  full  illumina- 
tion is  termed  the  effective  aperture.   Its  effect  will  be 
discussed  in  a  later  section.   Figure  11  shows  for  comparison 
the  image  plane  intensity  profile  of  a  fully  illuminating 
plane  wave. 

2.   Aperture  Size 

The  universal  feature  of  realistic  weighting  functions 
is  that  of  a  finite  aperture  corresponding  to  a  finite  period 
of  signal  x  represented  by  the  unity  amplitude  window  function 
[Ref.  5]. 

w(t,i)  =  rect  (|  -  |  )  (20) 

where 

T  =  D/Vs 

Early  approaches  to  acoustooptical  spectrum  analysis  modeling 
considered  the  input  signal  as  a  continuous  wave.   The  effec- 
tive aperture  of  the  selected  Bragg  deflector  truncated  the 
incident  light,  such  that  only  the  segment  of  acoustic  wave 
within  the  optical  aperture  region  was  intercepted.   That 
segment  of  acoustic  wave  corresponds  to  the  finite  period  of 
signal  being  analyzed. 

Improvements  in  Bragg  deflector  manufacturing  tech- 
niques have  yielded  high  quality,  low  acoustic  velocity 
deflector  crystals.   These  crystals  are  of  increasingly  larger 
dimensions.   The  Rayleigh  limit  resolution,  N,  is  defined  as 
the  total  deflection  angle  divided  by  the  angular  spot  size. 
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Figure   11. 


Image  Plane  Intensity  Profile  (Magnitude)  Produced 
by  a  Monotonic  7  0  MHz  Signal  with  an  Incident  Plane 
Wave  and  a  Simple  Transducer 
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The  relation  simplifies  to 

N  =  ti  x  if  (21) 

where 

ti  =  acoustic  transit  time  across  the  optical  aperture 

Af  =  deflector  bandwidth 
For  the  AOD-70  deflector  (assuming  that  full  aperture  illumi- 
nation is  provided) 

ti  =  2.56  x  D  (microsec)  (22) 

where  D  is  the  optical  aperture  in  centimeters. 
Using  the  AOD-70  parameters  of  a  40  MHz  bandwidth  and  a  maxi- 
mum tA  of  10  microsec  based  on  the  full  available  aperture 
of  3.9  cm,  calculation  yields  a  Rayleigh  resolution  capacity 
of  400  spots.   To  fill  the  full  aperture  anamorphic  optics 
such  as  cylindrical  lenses  or  telescopes  must  be  used.   Using 
standard  procedure  for  the  detector  array  utilization  each 
spot  is  focused  to  fall  across  two  adjacent  photosites.   The 
Rayleigh  resolution  capacity  of  400  spots  provides  for  a 
maximum  frequency  resolution  of  .1  MHz  for  the  full  40  MHz 
bandwidth. 

Using  a  Gaussian  beam  with  spot  size  of  1.945  milli- 
meters to  fully  suppress  the  sidelobes  in  the  image  plane 
intensity  profile  results  in  an  additional  penalty  in  a 
decrease  in  the  Rayleigh  resolution  capacity  of  the  "time- 
bandwidth  product,"  as  given  by  Equation  21.   The  interaction 
length  is  no  longer  the  Bragg  cell  aperture  but  rather  the 
effective  aperture,  in  this  case  estimated  as  the  region 
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twice  the  spot  size  symmetrical  about  the  cell  spatial  midpoint, 
This  results  in  a  smaller  optical  aperture,  a  smaller  acoustic 
transit  time,  and  a  resultant  decrease  in  the  Rayleigh  resolu- 
tion.  Thus,  while  a  choice  of  spot  size  of  this  magnitude 
will  suppress  sidelobes,  it  requires  a  possibly  unacceptable 
tradeoff  in  frequency  resolution. 

The  larger  apertures  bring  an  additional  consideration 
to  the  AOSA  model  construction.   A  typical  application  input 
signal  pulse  might  be  of  one  microsecond  in  duration  while 
the  total  acoustic  transit  time  for  the  AOD-7  0  is  ten  micro- 
seconds.  A  one  microsecond  pulse  input  to  the  deflector  will 
persist  in  terms  of  the  frequency  spectra  it  contains  and 
therefore,  in  the  deflected  beam  components  incident  on  the 
image  plane  for  the  full  ten  microsecond  travel  time.   This 
effect  of  "time  spreading"  of  the  r-f  signal  is  significant 
for  its  enhancement  of  the  detection  probability  for  short 
duration  signals.   In  the  application  of  AOSA  techniques  to 
Electronic  Support  Measures  (ESM)  receivers,  this  is  of  prime 
import.   A  0.1  microsecond  r-f  signal,  which  is  currently 
beyond  the  capability  of  many  traditional  technology  base 
receivers,  would  be  stretched  to  a  1.0  microsecond  pulse 
greatly  improving  its  detection  probability.   It  should  be 
noted  that  while  this  "time  spreading"  increases  detection 
probability,  it  reduces  the  capability  for  accurate  direct 
pulse  width  measurement,  also  a  highly  desirable  factor  in 
ESM  application. 
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In  summary  the  aperture  size  weighting  function 
imposes  three  constraints  on  the  AOSA  model.   The  physical 
parameter  of  the  optical  aperture  is  modeled  by  a  simple 
unity  amplitude  window  function.   The  Rayleigh  resolution 
and  thus  frequency  resolution  is  simulated  by  the  technique 
of  zero  padding  while  applying  the  Fourier  transforms  to 
achieve  a  digital  record  length  which  will  yield  the   appro- 
priate frequency  resoltuion.   Zero  padding  is  a  digital 
implementation  technique  used  in  Fourier  analysis  whereby 
the  desired  sampled  signal  data  is  entered  followed  by  a 
data  string  of  zeroes.   This  increases  the  record  length  to 
provide  a  desired  resolution  in  the  transformed  domain  with- 
out affecting  the  sampled  data  validity.   The  transformed 
domain  resolution  is  simply  the  inverse  of  the  record  length. 
The  finite  signal  length  is  modeled  with  a  unity  amplitude 
window  function  implicit  in  the  input  signal  as  it  is  digitally 
implemented.   Since  the  period  of  signal  analyzed  may  be 
shorter  than  the  acoustic  travel  time  of  the  optical  aperture, 
it  is  desirable  to  analyze  the  deflection  caused  by  inter- 
action of  the  pulse  acoustic  segment  at  various  points  as  it 
traverses  the  optical  aperture.   A  shiftable  digital  array 
was  utilized  controlled  by  a  time  delay  variable;  zero  delay 
corresponding  to  a  symmetric  straddling  of  the  acoustic  pulse 
segment  about  the  spatial  midpoint  of  the  optical  aperture. 
This  time  delay  variable  forms  the  basis  for  a  dynamic  AOSA 
program  to  be  discussed  in  a  later  section. 
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3.   Acoustic  Attenuation 

As  the  acoustic  wave  travels  the  length  of  the  Bragg 
deflector,  it  is  subject  to  acoustic  attenuation.   The  ampli- 
tude of  the  acoustic  wave  decreases,  with  a  resultant  decrease 
in  the  ability  to  interact  with  the  incident  light  wave. 
Hecht  [Ref.  5]  treats  the  acoustic  attenuation  as  a  spatial 
weighting  function  of  the  form; 

w(t,a,x)  =  exp(-at)  (23) 

where  a  is  a  frequency  dependent  loss  factor 

,-m 
a   f 

a   =   57T-T2 (24) 

20    log1Q   e 

with 

dB 


a      in 


ysec   GH 


TIT 


Values  for  a     have  been  tabulated  for  common  acoustooptic 
o 

modes.   Acoustic  attenuation  typically  follows  power  law 
frequency  dependence  with  exponent,  m,  near  2. 

The  AOD-70  Bragg  deflector  is  manufactured  with  high 
quality  flint  glass  as  the  interaction  medium  and  Indium 
bonded  Lithium  Niobate  for  the  piezo-electric  transducers. 

The  a  values  for  this  type  of  deflector  range  from  .1  to  .15 

o 
(dB/microsecond  GHZ  ) .   A  plot  of  the  more  stringent  condi- 
tion is  shown  as  Figure  12.   Calculation  of  the  acoustic 
attenuation  for  the  AOD-70  device  parameters  in  terms  of 
a  normalized  transmission  factor  was  carried  out.   The 
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results  are  plotted  as  Figure  13.   These  results  show  a  near 
unity  transmission  for  the  region  of  concern,  from  50  to 
90  MHz.   Since  the  transmission  loss  was  negligible  for  this 
application,  an  acoustic  attenuation  factor  of  Hecht ' s  form 
was  not  implemented  in  the  AOSA  programs. 

If  the  use  of  the  AOSA  programs  is  envisioned  in 
higher  frequency  ranges,  acoustic  attenuation  will  become  a 
significant  factor.   Figure  13  shows  increased  loss  as  the 
frequency  nears  1  GHZ.   Further  discussion  of  the  effect  of 
acoustic  attenuation  on  weak  truncations  (Gaussian  beam)  and 
strong  truncations  (uniform  optical  beam)  is  given  by  Hecht 
[Ref.  21]. 

4 .   Diffraction  Efficiency 

For  a  simple  transducer  the  Bragg  angle  cannot  be 
exactly  matched  by  the  incident  collimated  light  beam  for 
all  frequencies  over  the  bandwidth  of  concern.   As  a  result, 
for  Bragg  regime  deflection  the  intensity  of  the  diffracted 
beam  is  a  function  of  the  applied  frequency.   The  ratio  of 
the  diffracted  beam  intensity  to  the  intensity  of  the  incident 
beam  is  defined  as  the  diffraction  efficiency.   An  analysis 
of  Bragg  regime  deflection  was  conducted  by  Reagan  in  Ref.  2 
and  is  summarized  here. 

Hecht  [Ref.  5]  expresses  diffraction  efficiency  in 
a  simplified  form  as: 


1  /     ~1 ~ 2 


*i       2  /-»  At  sin' 
— -  =   sin   (Tl  '  )  


1 1  (FFm"F2) 


SS  (FF  -F2) 
2  L    m 
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with 


and 


2  6     2  PL 

q   =     JLln!£l_^ (26) 

2A2      pv3     4    cos2    9 
o  o 


L      =   nA2    cos    0   /A  (27) 

o  o  o'    o 


where:   Ii  =  diffracted  beam  intensity 

I   =  undeflected  beam  intensity 

o  x 

A  =  free  space  laser  wavelength 

n  =  refractive  index 

p  =  elastooptic  coefficient 

p  =  mass  density 

A  =  acoustic  wavelength 

v  =  acoustic  wave  velocity 

P  =  acoustic  power 
a 

H  =  acoustic  beam  height 

9   =  Bragg  deflection  angle 

L  =  interaction  length 

L  =  characteristic  interaction  length 

F  =   normalized  center  frequency  =  1 

F   =  input  frequency  value  normalized  to  the  center 
frequency  value  f 

Equation  26  shows  the  quantity  r\    to  be  constant  with 
respect  to  the  input  signal  frequency.   Thus  the  initial  sin 
squared  term  of  Equation  25  may  be  ignored  when  considering 
the  bandshape  of  a  normal  acoustooptic  beam  deflector.   A  plot 
of  Equation  25  for  various  L/L   ratios  is  shown  as  Figure  14. 
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Figure  14.   Diffraction  Efficiency  Bandshape  (after  Ref.  2) 
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The  3dB  bandwidth  is  given  by  Equation  28  and  is  also  from 

Ref.  5. 

Af  =  1.8  f   L  /L  (28) 

o   o 

For  application  of  a  Bragg  regime  deflector  as  an 
acoustooptical  spectrum  analyzer,  it  is  desirable  to  maintain 
a  bandshape  as  near  unity  and  as  uniform  as  possible  over 
the  bandwidth  of  concern.   The  characteristic  interaction 
length  is  shown  by  Equation  27  to  be  constant  with  respect  to 
dimensions  of  the  Bragg  deflector.   The  trend  for  AOSA 
applied  Bragg  deflectors  is  that  of  increased  size  with  an 
attendant  increased  interaction  length  L.   This  reflects 
both  improvement  in  manufacturing  techniques  permitting  the 
production  of  high  quality  large  deflectors  and  the  desire 
for  increased  resolution  made  possible  with  larger  inter- 
action length.   The  increased  interaction  length,  however, 
results  in  higher  value  L/L,   ratios  with  the  observed 
effect  of  narrowed  bandwidths. 

Attempting  to  increase  bandwidth  by  shortening  the 
interaction  length  is  of  limited  benefit  since  the  diffraction 
shifts  from  the  Bragg  regime  to  the  Raman-Nath  regime. 
Raman-Nath  regime  diffraction  is  unsuitable  for  this  appli- 
cation as  it  supports  many  higher  order  diffraction  modes 
resulting  in  interference  created  by  additional  beams. 
Increasing  the  refractive  index  or  elastooptic  coefficient 
of  the  interaction  medium,  decreasing  the  acoustic  beam 
height,  or  decreasing  the  acoustic  velocity  or  the  mass 
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density  of  medium  can  improve  the  diffraction  efficiency  to 
a  limited  extent.   An  increase  of  acoustic  beam  power  will 
increase  diffraction  efficiency.   The  increased  power  also 
tends  to  heat  the  interaction  medium,  creating  temperature 
gradients  which  degrade  resolution  and  information  capacity. 
Increased  difficulty  in  internal  reflection  damping  is 
experienced,  permitting  unwanted  signal  persistence. 

The  only  effective  manner  in  which  to  increase  the 
bandwidth  is  to  increase  the  cell's  center  frequency  f  . 
The  maximum  frequency  which  the  input  transducer  can  effec- 
tively couple  to  the  medium  is  determined  primarily  by  the 
bonding  technique  used  to  attach  the  piezo-electric  strips 
to  the  interaction  medium.   Thus  the  ability  to  increase  the 
center  frequency  is  dependent  on  the  input  signal  transducer. 
Normally  Bragg  deflectors  are  limited  to  operation  over  one 
octave . 

In  summary,  for  a  simple  transducer  little  can  be 
done  to  improve  diffraction  efficiency  or  increase  bandwidth. 
Beam  steering,  however,  can  significantly  improve  both 
diffraction  efficiency  and  bandwidth.   Commonly  employed 
beam  steering  techniques  include  first  order  planar  and 
stepped  arrays.   Figure  15  shows  a  typical  first  order 
stepped  array.   Pinnow  [Ref.  22]  discussed  the  relative 
merits  of  simple,  planar  and  stepped  array  transducers. 
The  larger  dimensioned  Bragg  cells  deflectors  are  almost 
always  of  the  stepped  array  type.   The  AOD-70  deflector  is 
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Figure  15.   First-order  Stepped  Array  Transducer 
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a  five  element  stepped  array.   Figure  16  is  the  manufacturer's 
supplied  efficiency  versus  frequency  response  curve. 

The  AOSA  program  has  both  a  simple  transducer  and 
a  stepped  array  approximation  option.   The  simple  transducer 
response  is  computed  in  accordance  with  isotropic  theory 
based  on  the  AOD-70  device  full  dimensions.   This  permits  a 
comparison  with  the  actual  stepped  array  response.   Figure  17 
shows  the  computed  diffraction  efficiency  curves  for  various 
L/L   ratios.   These  were  obtained  by  sweeping  a  single  varia- 
ble frequency  over  the  band  of  concern  and  extracting  the 
magnitude  value  for  the  instantaneous  frequency  applied.   In 
comparison  with  Figure  14,  which  contains  the  bandshapes  for 
the  diffraction  efficiency  theory ,  a  rippling  is  noticeable 
at  the  lower  frequencies.   This  rippling  is  an  undesirable 
artifact  of  the  digital  Fourier  transform  application.   The 
ripple  amplitude  is  comparatively  small  and  falls  outside  of 
region  of  concern;  hence  it  was  not  corrected. 

The  stepped  array  bandshape  that  was  included  in  the 
model  is  an  approximation  based  on  the  near  uniform  response 
that  was  experimentally  measured  for  the  AOD-70  device.   A 
slight  roll-off  in  response  occurs  at  the  endpoints  of  the 
band,  but  for  a  first  order  approximation  a  uniform  window 
function  utilizing  the  mean  efficiency  was  employed. 

To  compare  the  response  of  the  models  of  the  two 
transducers,  an  input  signal  consisting  of  an  in-phase 
frequency  comb  with  5  MHz  spacing  was  applied.   Figure  18  shows 
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Figure  16.   AOD-70  Bragg  Deflector  Efficiency  Curve 
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Figure  17.   AOSA  Implemented  Diffraction  Efficiency  Bandshape 
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Figure    18. 


Image  Plane  Intensity  Profile  (Magnitude)  Produced 
by  a  5  MHz  Incremental  Comb  Signal  with  an  Incident 
Plane  Wave  and  a  Simple  Transducer 
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the  simple  transducer  response.   Figure  19  shows  the  stepped 
array  response.   Comparison  shows  severe  attenuation  of  the 
frequency  components  beyond  80  MHz  and  a  lesser  degree  of 
attenuation  from  50  to  60  MHz  for  the  simple  transducer. 
On  the  other  hand  the  stepped  array  transducer  shows  a  very 
uniform  response  but  with  a  tradeoff  in  decreased  overall 
intensity  magnitude. 

5.   Transducer  Efficiency 

The  measure  of  the  transducer's  ability  to  convert 
the  r-f  input  signal  components  to  representative  acoustic 
waves  is  termed  the  transducer  efficiency.   The  transducer 
efficiency  is  dependent  upon  the  transducer  material,  the 
bonding  technique  used  to  fasten  the  transducer  to  the  inter- 
action medium,  and  the  transducer  impedance  matching.  A  funda- 
mental tradeoff  relationship  exists  between  the  flatness  of 
the  frequency  response  and  the  amount  of  power  delivered  by 
the  transducer  within  the  desired  band.   The  power  levels 
used  to  drive  the  transducer  are  relatively  low,  thus  we  can 
afford  a  power  conversion  inefficiency  and  seek  as  uniform 
a  response  as  possible.   All  transducers  operate  at  less  than 
unity  efficiency;  however,  for  purposes  of  the  model,  it  was 
assumed  that  the  desired  uniform  response  was  of  paramount 
importance.   An  ideal  case  would  provide  a  transducer  efficiency 
response  with  little  power  attenuation  and  a  uniform  amplitude 
for  all  frequencies  within  the  range  of  the  system  application. 
Current  technology  has  achieved  near  uniform  response  bandshapes 
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Figure  19.   Image  Plane  Intensity  Profile  (Magnitude)  Produced 
by  a  5  MHz  Incremental  Comb  Signal  with  an  Incident 
Plane  Wave  and  a  Stepped  Array  Transducer 
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over  bandwidths  of  70  MHz,  as  is  present  in  the  AOD-70  cell. 
Chang  [Ref.  23]  discusses  the  state  of  current  transducer 
technology. 

Based  on  the  above  the  AOSA  model  was  implemented 
with  a  unity  amplitude  window  for  the  transducer  efficiency 
function.   The  program  structure  permits  easy  input  of  non- 
unity  bandshapes  should  this  be  desired.   Extreme  broadband 
applications  of  acoustooptical  spectrum  analysis  will  result 
in  nonuniform  frequency  responses  requiring  such  modifications 
before  using  the  AOSA  program. 

D.   AOSA  MODEL  UTILIZATION  AND  EXTENSIONS 

The  AOSA  program  was  designed  to  model  direct  detection 
acoustooptical  spectrum  analyzer  systems  in  general.   In  its 
current  configuration  it  utilizes  the  component  parameters 
of  the  NPS  laboratory  system  with  the  addition  of  the  AOD-70 
Bragg  cell.   The  coding  for  the  Fortran  program  was  selected 
to  maintain  as  much  generality  as  possible,  so  as  to  permit 
application  for  analysis  under  other  choices  of  components. 
An  attempt  was  made  to  minimize  the  number  of  preliminary 
calculations  a  prospective  user  would  have  to  perform  prior 
to  utilization  of  the  program.   For  example,  the  program 
computes  the  Fourier  transform  lens  focal  distance  given  the 
deflector  physical  parameters,  the  laser  characteristics, 
and  the  desired  detector  array  spot  size.   A  list  of  the 
required  input  physical  parameters  is  given  in  Table  II. 
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TABLE  II 
AOSA  Modeling  Parameters 

DCELL  =  effective  optical  aperture  of  Bragg  cell 

LW     =  Bragg  cell  interaction  length 

VS     =  acoustic  velocity  of  interaction  medium 

REFIND  =  refractive  index  of  interaction  medium 

LAMBDA  =  laser  wavelength 

BW     =  beam  width  of  collimated  Gaussian  beam  (3  dB)  at 

the  Bragg  cell  aperture 
SPOT   =  desired  detector  array  spot  size 
PULSE  =  r-f  signal  pulse  duration 
TDELAY  =  time  delay  for  examining  signal  pulse  at  varied 

points  during  the  acoustic  transit  period 
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The  remaining  program  variables  are  primarily  concerned  with 
control  of  the  program  including  output  formats  and  the  var- 
ious subroutines  utilized  within  the  program. 

The  main  exception  to  the  above  is  the  created  variable 
TDELAY,  for  time  delay,  which  permits  the  user  to  examine 
the  representative  acoustic  pulse  at  frozen  intervals  of  time 
during  its  transit  of  the  optical  interaction  aperture.   This 
suggested  the  possibility  of  adding  another  dimension  to  the 
AOSA  program,  to  create  in  effect  a  dynamic  AOSA  model  which 
would  examine  the  response  of  the  acoustooptic  spectrum 
analyzer  periodically  by  looping  through  the  basic  static 
algorithm  with  decreasing  time  delays.   The  results  are 
presented  in  a  three-dimensional  isometric  or  perspective 
projection.   This  permits  the  user  to  simulate  the  response 
of  the  image  plane  intensity  as  the  acoustic  pulse  transits 
segments  of  the  optical  interaction  region  rather  than  at  a 
single  fixed  point  of  time.   To  simplify  the  graphic  output 
the  present  configuration  of  the  dynamic  AOSA  program  (called 
TAOSA)  is  set  to  examine  the  first  fraction  of  a  microsecond 
as  the  acoustic  pulse  enters  the  optical  region.   Computer 
simulations  were  conducted  over  the  entire  transit  period 
and  to  verify  the  symmetry  of  response  as  the  acoustic  pulse 
exited  the  optical  region. 

Both  the  static  and  dynamic  programs  provide  the  user 
with  options  to  simulate  various  system  configurations. 
Types  of  transducers,  incident  beam  profiles,  and  output  in 
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either  relative  magnitude  or  decibel  form  is  available.   The 
static  program  offers  a  choice  of  several  levels  of  graphical 
output,  from  a  coarse  line  plot  to  two  Versatec  offline 
plotting  options.   Currently  both  programs  utilize  a  deter- 
ministic form  for  the  input  r-f  signal.   The  programs  could 
be  easily  modified  to  accept  sampled  r-f  data. 

The  following  series  of  figures  are  typical  of  the  outputs 
of  the  AOSA  and  TAOSA  programs.   They  were  selected  to  demon- 
strate the  ability  of  the  programs  to  model  a  variety  of 
complexity  in  input  signals,  to  demonstrate  the  quality  of 
resolution  achievable  by  an  acoustooptical  spectrum  analyzer 
and  for  comparison  with  the  known  spectral  content  of  the 
inputed  signals. 

Figure  20  is  the  input  r-f  signal  magnitude  for  a  70  MHz 
monotone.   Figures  21  and  22  are  outputs  using  the  decibel 
option  of  the  AOSA  program.   Figure  21  shows  the  relative 
sidelobe  levels  by  interaction  of  the  monotone  with  an  inci- 
dent plane  wave.   For  comparison  Figure  22  shows  the  decrease 
in  the  relative  sidelobe  level  produced  by  the  monotone  with 
an  incident  Gaussian  beam  of  1.045  millimeter  spot  size. 
Figures  23  and  24  show  the  dynamic  response  of  the  same 
monotone  with  the  two  wave  types.   The  x  axis  on  all  TAOSA 
program  graphics  represents  time  as  the  acoustic  pulse  transits 
the  Bragg  cell  aperture  region.   Time  zero  corresponds  to  the 
instant  when  the  leading  edge  of  the  acoustic  wavetrain 
crosses  into  the  aperture  region.   Figure  2  3  shows  the  dynamic 
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Figure  20.   Input  RF  Signal  Magnitude  for  a  70  MHz  Monotone 
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Figure  21.   Image  Plane  Intensity  Profile  (Decibel)  Produced 
by  a  Monotonic  7  0  MHz  Signal  with  an  Incident 
Gaussian  Beam  and  a  Stepped  Array  Transducer 
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Figure  22.   Image  Plane  Intensity  Profile  (Decibel)  Produced 
by  a  Monotonic  7  0  MHz  Signal  with  an  Incident 
Plane  Wave  and  a  Stepped  Array  Transducer 
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Figure  23. 


Image  Plane  Intensity  Profile  (Magnitude)  Produced 
by  a  Monotonic  70  MHz  Signal  with  an  Incident  Plane 
Wave  and  a  Stepped  Array  Transducer  (Dynamic  Model) 
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response  of  the  monotone  with  an  incident  plane  wave. 
A  much  slower  response  is  seen  in  Figure  24  showing  the 
response  produced  by  an  incident  Gaussian  wave  of  1.945 
millimeter  spot  size.   The  Gaussian  response  is  shown  to 
occur  later  in  the  cell  aperture  transit  period  and  to  last 
for  a  relatively  short  time  when  compared  to  the  plane  wave 
response.   The  symmetric  roll-off  is  already  visible  at  the 
rightmost  edge  of  the  figure. 

Figure  25  is  the  input  r-f  signal  magnitude  for  a  simple 
sinusoidally  modulated  FM  signal  with  6=8.0  and  f   =  5  MHz. 
6  is  the  modulation  index  and  is  defined  as 

B«r  (29) 

m 

Figure  26  shows  the  image  plane  intensity  profile  produced 
by  the  interaction  of  the  signal  shown  in  Figure  25  with  an 
incident  plane  wave.   Figure  27  shows  the  image  plane  inten- 
sity profile  produced  by  the  same  signal  interacting  with  an 
incident  Gaussian  beam  of  1.94  5  millimeter  spot  size. 
Figures  28  and  29  are  the  analagous  dynamic  program  outputs. 
Again  the  Gaussian  beam  interaction  shows  a  slower  and 
smaller  response.   The  advantage  of  sidelobe  suppression, 
however,  is  apparent  in  the  extremely  clear  presentation  of 
the  FM  intermodulation  products. 

Figure  30  shows  the  image  plane  intensity  profile  produced 
by  a  linear  chirp  with  an  incident  plane  wave.   Figure  31 
shows  the  profile  produced  by  the  interaction  with  a  Gaussian 
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Figure  24. 


Image  Plane  Intensity  Profile  (Magnitude)  Produced 
by  a  Monotonic  70  MHz  Signal  with  an  Incident 
Gaussian  Beam  and  a  Stepped  Array  Transducer 
(Dynamic  Model) 
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Figure  25.   Input  RF  Magnitude  for  an  FM  Modulation  Signal 

with  8  =  3.0,  f   =  5  MKz 
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Figure  26.   Image  Plane  Intensity  Profile  (Magnitude)  Produced 
by  an  FM  Modulation  Signal,  6  =  8.0,  f   =5  MHz, 
with  an  Incident  Plane  Wave  and  a  Stepped  Array 
Transducer 
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Figure    27. 


Image  Plane  Intensity  Profile  (Magnitude)  Produced 
by  an  FM  Modulation  Signal,  3  =  8.0,  f   =5  MHz, 
with  an  Incident  Gaussian  Beam  and  a  Stepped  Array 
Transducer 
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Figure  28. 


Image  Plane  Intensity  Profile  (Magnitude)  Produced 
by  an  FM  Modulation  Signal,  3  =  8.0,  f   =  5  MHz, 
with  an  Incident  Plane  Wave  and  a  Stepped  Array 
Transducer  (Dynamic  Model) 
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Figure  29. 


Image  Plane  Intensity  Profile  (Magnitude)  Produced 
by  an  FM  Modulation  Signal,  3=  8.0,  f   =5  MHz, 
with  an  Incident  Gaussian  Beam  and  a  Stepped  Array 
Transducer  (Dynamic  Model) 
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Figure    30. 


Image  Plane  Intensity  Profile  (Magnitude)  Produced 
by  a  Linear  Chirp  with  Af   =80  MHz,  f  =50  MHz, 
t  =  2  ysec  with  an  Incident  Plane  Wave°and  a 
Stepped  Array  Transducer 
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Figure  31.   Image  Plane  Intensity  Profile  (Magnitude)  Produced 
by  a  Linear  Chirp  with  Af   =  80  MHz ,  f   =  50  MHz, 
t  =  2  ysec  with  an  Incident  Gaussian  B§am  and  a 
Stepped  Array  Transducer 
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beam.   For  both  the  chirp  parameters  are  Af   =80  MHz, 

m 

f   =50  MHz.  and  t  =  2  microseconds.   Both  show  high  and 

o  3 

low  frequency  ripples.   The  high  frequency  ripple  is  an  arti- 
fact of  the  digital  implementation  of  the  linear  chirp.   The 
lower  frequency  ripple  present  is  the  effect  of  the  funda- 
mental and  sidelobe  interaction.   The  Gaussian  interaction 
shows  the  effects  of  its  inherent  slower  response  and  degraded 
frequency  resolution  in  the  attenuation  of  magnitude  over  the 
lower  half  band  frequency  components. 

Figure  32  shows  the  dynamic  response  of  the  same  linear 
chirp  and  an  incident  plane  wave.   The  severity  of  the  inter- 
action of  the  fundamental  and  sidelobe  components,  as  addi- 
tional frequency  components  enter  the  Bragg  cell,  manifests 
itself  in  a  dynamic  transient  response.   Only  after  a  con- 
siderable interval  does  the  response  settle  into  a  steady 
state  displaying  the  expected  chirp  spectral  profile.   For 
comparison  Figure  33  shows  the  same  linear  chirp  interacting 
with  an  incident  Gaussian  beam  of  1.94  5  millimeter  spot  size. 
In  this  case  the  sidelobe  suppression  is  a  distinct  advantage 
producing  a  smooth  steady  state  response  almost  from  the 
onset  of  the  interaction.   It  must  be  kept  in  mind  that  the 
interaction  will  continue  to  begin  much  later  than  that  of 
the  plane  wave  and  thus  the  plane  wave  will  achieve  steady 
state  earlier  than  the  Gaussian. 

The  TAOSA  program  also  has  the  decibel  output  option. 

Figure  34  shows  the  relative  decibel  image  plane  intensity 

of  a  monotone  of  70  MHz  interacting  with  a  plane  wave. 
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Figure  32. 


Image  Plane  Intensity  Profile  (Magnitude)  Produced 
by  a  Linear  Chirp  with  Af   =80  MHz,  f   =50  MHz, 
t  =  2  ysec,  with  an  Incident  Plane  Wave  and  a 
Stepped  Array  Transducer  (Dynamic  Model) 
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Figure  33.   Image  Plane  Intensity  Profile  (Magnitude)  Produced 

by  a  Linear  Chirp  with  Af   =  80  MHz,  f   =  50  MHz, 

t  =  2  usee  with  an  Incident  Gaussian  Beam  and  a 
Stepped  Array  Transducer  (Dynamic  Model) 
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Figure  34. 


Image  Plane  Intensity  Profile  (Decibel)  Produced 
by  a  Monotonic  70  MHz  Signal  with  an  Incident 
Plane  Wave  and  a  Stepped  Array  Transducer  (Dynamic 
Model) 
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III.   SUMMARY  AND  CONCLUSIONS 

A.  SUMMARY 

This  thesis  has  presented  an  analysis  of  Bragg  diffrac- 
tion theory  and  direct  detection  acoustooptical  spectrum 
analysis  modeling  techniques.   Based  on  the  study  of  the 
various  modeling  theories  a  suitable  approach  for  a  first 
order  simulation  was  selected  and  used  to  develop  a  Fortran 
program  modeling  the  direct  detection  acoustooptical  spectrum 
analysis  process.   Extensions  of  the  basic  theory  have  been 
added  to  include  the  factors  of  diffraction  efficiency, 
transducer  efficiency,  irradiance  profiles  of  the  incident 
laser  illumination,  the  aperture  size  of  the  Bragg  cell, 
and  acoustic  attenuation  of  the  generated  acoustic  wave. 
A  second  program  dynamically  models  the  AOSA  process  incor- 
porating a  time  variable  to  simulate  the  response  of  the  AOSA 
as  the  acoustic  wave  transits  the  interaction  region.   The 
programs  were  used  to  model  the  NPS  laboratory  AOSA  using 
the  parameters  from  that  system. 

B.  AREAS  FOR  FURTHER  STUDY 

Further  effort  related  to  this  thesis  is  required  in  the 
following  areas. 

1.   Nonlinearities 

The  theoretical  basis  for  the  AOSA  model  program  does 
not  take  into  account  the  nonlinear  effects  of  crossmodulation 

80 


(the  compression  of  one  signal  due  to  another)  or  inter- 
modulation  (the  interaction  of  two  or  more  signals) .   While 
the  restriction  of  the  operation  of  the  AOSA  to  the  Bragg 
regime  eliminates  much  of  the  effects  of  these  nonlineari- 
ties,  the  overall  effect  on  the  replication  efficiency  is 
still  at  issue.   Further  analysis  is  needed  to  determine 
if  these  factors  seriously  degrade  the  model's  accuracy. 
Hecht  discusses  this  issue  in  Ref.  21. 

2.  Time  Varying  Signals 

The  ability  of  the  AOSA  to  respond  to  rapid  time 
varying  signals  (which  will  depend  on  the  speed  of  response 
of  the  detector  array  and  the  Bragg  cell)  needs  to  be 
examined.   Currently  the  TAOSA  program  provides  for  a  fixed 
delay  for  the  simulation  of  overall  system  response  time. 

3 .  System  Components 

The  physical  response  and  performance  of  the  AOSA 
systems  components  warrant  further  study.   Currently  only 
a  gross  simulation  of  the  system  response  time  can  be  made 
by  insertion  of  a  fixed  time  delay.   Closer  examination  of 
the  speed  of  the  detector  array  with  particular  emphasis  on 
the  uniformity  of  photosite  response  is  needed.   Integrated 
optical  devices  for  AOSA  and  AOSE  applications  will  use 
semi-conductor  lasers.   These  lasers  will  need  to  be 
examined  for  the  system  effects  produced  by  fluctuations  in 
intensity  and  variations  in  coherence.   The  current  HeNe 
laser  used  in  the  NPS  system  should  be  examined  for  its 
intensity  and  coherence  fluctuations  and  resultant  effects. 
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C.   CONCLUSIONS 

Comparison  of  the  AOSA  model  image  plane  intensity 
profiles  with  known  AOSA  system  output  spectral  profiles 
demonstrates  the  AOSA  program  as  a  satisfactory  first  order 
simulation.   The  acoustooptical  spectral  analysis  process 
yields  high  quality  replication  of  the  spectral  content  on 
input  signals.   Future  coupling  of  the  AOSA  program  with 
optical  excision  routines  previously  developed  may  provide 
additional  insight  into  the  application  of  acoustooptical 
spectral  analysis  and  excision  for  use  in  wideband,  high 
density  signal  environments. 
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APPENDIX  A 
ACOUSTOOPTIC  SPECTRUM  ANALYZER  (AOSA)  PROGRAM 

This  Fortran  program  calculates  the  relative  Bragg 
deflected  image  plane  intensity  profiles  for  multifrequency 
input  r-f  signals  to  a  modeled  acoustooptic  spectrum  analyzer 
(AOSA) .   The  physical  parameters  used  in  the  modeling  of  the 
AOSA  are  discussed  within  the  program  and  are  summarized  in 
Table  II  of  the  thesis. 

The  program  provides  several  additional  selectable 
variables : 

ITEST   -  provides  for  graphical  and  numerical  or 

graphical  only  output 
IWAVE   -  provides  for  choice  of  plane  or  Gaussian 

incident  light  waves 
ITRANS  -  provides  for  choice  of  simple  or  first 

order  stepped  array  transducer  simulations 
IDB    -  provides  for  choice  of  graphical  output  in 

magnitude  or  relative  dB  formats 
IGRAF   -  provides  for  choice  of  line  printer  or 

Versatec  graphical  plotting 
IPLOT   -  provides  for  choice  of  PLOTG  or  DRAWP 
Versatec  plotting  subroutines 
Use  of  these  variables  is  discussed  within  the  program. 
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Additional  alterable  variables  exist  within  the  program 
related  to  the  Fourier  transform  process.   It  is  believed 
that  the  present  configuration  will  suffice  for  most  general 
applications. 

This  program  utilizes  the  PLOTG  and  DRAWP  Versatec 
software  subroutines,  the  line  plot  subroutine  PLOTP,  and 
FOURT,  which  calculates  the  forward  and  inverse  Cooley-Tukey 
Fast  Fourier  Transform  of  multidimensional  complex  data. 
Prospective  utilizers  at  the  Naval  Postgraduate  School  should 
refer  to  the  W.  R.  Church  Computer  Center  for  detailed 
information. 
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APPENDIX  B 
DYNAMIC  ACOUSTOOPTIC  SPECTRUM  ANALYZER  PROGRAM 

The  dynamic  acoustooptic  spectrum  analyzer  program  (TAOSA) 
is  essentially  a  variation  of  the  AOSA  Fortran  program.   It 
computes  the  time  delay  to  shift  the  acoustic  pulse  outside 
of  the  optical  interaction  aperture  and  then  incrementally 
decreases  this  delay  to  simulate  the  action  of  the  pulse  as 
it  transits  the  optical  interaction  region.   The  system 
response  time  may  be  simulated  by  adding  a  fixed  time  delay. 

The  variables  used  in  the  program  are  discussed  both 
within  the  program  and  in  Appendix  A.   An  additional  variable, 
TINT,  is  added  as  an  acceleration  factor  that  varies  the  in- 
crement size   based  on  the  level  of  intensity  of  the  incident 
laser  beam. 

The  results  are  presented  as  a  three-dimensional  isometric 
or  perspective  plot.   Perspective  users  at  the  Naval  Post- 
graduate School  should  refer  to  the  W.  R.  Church  Computer 
Center  for  detailed  information  on  the  three-dimensional 
plotting  routine  PLT3D1  in  use. 
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APPENDIX  C 

ARBITRARY  PROFILE  SOLUTION  DETAILS 
(after  Ref.  15) 

The  complex  transmission  coefficient  T(cJ>)  is  defined  as 

T(4>)  =  4k*   fa   -  a    )  (A.  m    -  B,  ;  ex?(-jk?L)/Det  (ci) 

Is      z  i    Da  Da  3x 

where 


Det  »  (a.   A    -aB   Jfa  A,  ,  -  a  B,,J 
2  aa     i  aa    i    dd  z    dd 


*^*(kab    '    *ab)(Xba   '    h 


)  (C2) 

a 


«i  =  Jjpr l/£,  -  #  +  [(/ft,  -  ^)2  +  (MkbW*),  (C3 ) 

a2  "  7717"  l^i  - ®  "  l<#i  -  ^)2  +  (MH^J2)V2\,  ( C4 ) 

x  [exp  (-cA^L)  -  R!fiJ|| exp(iASrL)j,  (C5) 

X  [exp<-iAk£)-rt$ffl$Jexp<iA  $,£>!.  (C6) 

8y*;"*['    p.(r-a.6.     ij-  1.2.3.  (C7) 

*  u  *  fcoV^coad  =»  V(2WX)iJ<1  -  3?  (C8) 

*t-(*l«i-i8ii)l/a  (C9) 

*a2,  - 1(1  -  (%)MatJ*Se,  - p^»«  (CIO) 

*L  -  1(1  -  (V2)iVfa?)^  -  $ll/2.  (Cll ) 

*f«  -  <*3«3  - /3D1/2  (C12) 

fcL-(*|ts-0ll)l/a  (CI  3) 

M  =  index  of  Modulation  of  the  permittivity  of  the  Bragg 
cell  medium 
z  ^relative  permittivity  of  air  *  1 
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£,  =  relative  permittivity  of  the  Bragg  cell  medium 


Z  £ 

3     a         1 

9      =   output  angle 
k        =*    2tt/X 

0 

3        =    C2tt/X)  /e"  sin    9  (C14) 

B        -      6      -    2TT/d  (C15) 

-l  0 
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